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Summary 

Coronaviruses are a diverse group of viruses that infect mammals and birds. Bats 
are reservoirs for several different coronaviruses in the Alphacoronavirus and 
Betacoronavirus genera. They also appear to be the natural reservoir for the ances¬ 
tral viruses that generated the severe acute respiratory syndrome coronavirus and 
Middle East respiratory syndrome coronavirus outbreaks. Here, we detected coro¬ 
navirus sequences in next-generation sequence data created from Eonycteris 
spelaea faeces and urine. We also screened by PCR urine samples, faecal samples 
and rectal swabs collected from six species of bats in Singapore between 2011 and 
2014, all of which were negative. The phylogenetic analysis indicates this novel 
strain is most closely related to lineage D Betacoronaviruses detected in a diverse 
range of bat species. This is the second time that coronaviruses have been detected 
in cave nectar bats, but the first coronavirus sequence data generated from this 
species. Bat species from which this group of coronaviruses has been detected are 
widely distributed across SE Asia, South Asia and Southern China. They overlap 
geographically, often share roosting sites and have been witnessed to forage on 
the same plant. The addition of sequence data from this group of viruses will 
allow us to better understand coronavirus evolution and host specificity. 


Introduction 

Over 70% of emerging or re-emerging infectious diseases 
originate in animals, with zoonotic RNA viruses responsi¬ 
ble for the majority of these cross-species spillover events 
(Woolhouse et al., 2005). Zoonotic outbreaks are amplified 
or dampened by a number of ecological, environmental 
and anthropogenic factors (Karesh et al., 2012). Change in 
land use is a principal driver for emerging infectious dis¬ 
eases by modifying the wildlife-human interface (Jones 
et al., 2008). Increase in human populations has caused 
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wild animals to lose their natural habitat and diminished 
resource availability, while bushmeat hunting and con¬ 
sumption have presented further opportunities for cross¬ 
species exposure (Wolfe et al., 2005). 

In South-East Asia, severe habitat loss has impacted the 
natural ecology of bats and occasionally resulted in zoono¬ 
tic virus spillover, most notably in the Nipah outbreaks in 
Malaysia and Bangladesh (Chua et al., 2000; Hsu et al., 
2004). Bats are reservoirs for many notable pathogenic 
viruses, including Hendra virus in Australia (Halpin et al., 
2000), the recent outbreaks of Ebola virus in sub-Saharan 
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Africa (Pourrut et al., 2009; Olival and Hayman, 2014; 
Ogawa et al., 2015), rabies and several other lyssaviruses 
(Banyard et al., 2014). Furthermore, zoonotic coron- 
aviruses (e.g. SARS-CoV and MERS-CoV) and the human 
coronavirus HCoV-229E, a causal agent of the common 
cold, appear to have resulted from zoonotic transfer from 
bats, either directly or through an intermediate amplifying 
host (Drexler et al., 2014). Indeed, bat coronaviruses are 
exceptionally diverse and genetic evidence provides further 
support that group 1 and group 2 mammalian coron¬ 
aviruses originated in bats (Vijaykrishna et al., 2007; Hu 
et al., 2015), highlighting the importance of bats as a 
potential source of new human diseases. 

Coronaviruses cause a wide variety of symptoms, with 
human coronaviruses causing respiratory infections in 
humans, while other animals suffer respiratory, gastroen¬ 
teric and more severe manifestations of disease (Peiris 
et al., 2003; Saif, 2004). Coronaviruses can be highly patho¬ 
genic, especially when they spill over into incidental hosts, 
as evidenced by the high fatality rates seen in the severe 
acute respiratory syndrome (SARS) in 2003 and the ongo¬ 
ing Middle East respiratory syndrome (MERS) outbreak, 
9.6% and 40%, respectively (Zumla et al., 2015). 

Coronaviruses are enveloped single-stranded, positive 
sense RNA viruses that belong to the Order Nidovirales , the 
Family Coronaviridae and the subfamily Coronavirinae 
(de Groot et al., 2012). There are four genera: Alphacoron- 
aviruSy Betacoronavirus , Deltacoronavirus and Gammacoron- 
avirus (Adams and Carstens, 2012). Viruses of the 
Deltacoronavirus and Gammacoronavirus genera primarily 
infect birds, with avian infectious bronchitis virus a promi¬ 
nent member of the latter that causes systemic infections of 
chickens, resulting in major economic losses (Jackwood, 
2012). Alphacoronavirues infect humans (HCoV-229E and 
HCoV-NL63), pigs, felines, canines and a number of different 
bat species (Drexler et al., 2014). An alphacoronavirus, that 
forms a sister group to HCoV-229E, has also been recently 
found to commonly infect camels (Sabir et al., 2016). Beta- 
coronaviruSy which is further classified into lineages A-D, is 
found in a wide variety of hosts including multiple bat 
species, humans and numerous domestic and peridomestic 
animals (Woo et al., 2009). Of note, this genus contains 
both SARS-CoVs (lineage B) and MERS-CoVs (lineage C). 

The large genomes and high rates of recombination wit¬ 
nessed in coronaviruses are suggested reasons for the ability 
of coronaviruses to jump into other species (Woo et al., 
2009). Interestingly, coronaviruses have some proof-read¬ 
ing mechanisms and are not subject to the infidelity of the 
RNA polymerase reading mechanism like other RNA 
viruses, such as influenza (Drake, 1993; Minskaia et al., 
2006). There is a deep phylogenetic divergence between the 
avian coronavirus-specific genera and the mammalian 
coronaviruses and a lack of evidence of avian-mammalian 
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cross-species transmission (Woo et al., 2012a); however, 
there is prolific evidence of these events occurring between 
different mammal species (Lau et al., 2012) and between 
different bird species (Chu et al., 2011). 

Detecting viruses that possess the capacity to effectively 
transmit across species greatly enhances our knowledge on 
the evolution and ecology of coronaviruses. Little is known 
about bat-virus interactions in Singapore, and although 
this country lacks endemic bat species, it shares the geo¬ 
graphical range of several South-East Asian bat species. 
Here, we use both RT-PCR and next-generation sequencing 
methods to detect coronaviruses in the bats of Singapore 
and examine their phylogenetic relatedness. 

Materials and Methods 

Sample collection 

Fresh faecal and urine samples were collected from a colony 
of Eonycteris spelaea on plastic sheets every 2 weeks from 
April 2011 until June 2015. Samples were placed into virus 
transport media at 4°C until transport to Duke-NUS where 
they were placed into a — 80°C freezer until screening. Addi¬ 
tionally, a total of 431 individual bats from six species: 
E. spelaea ( n = 169), Cynopterus brachyotis (n = 144), Pen- 
thetor lucasi ( n = 79), Macroglossus minimus (n = 2), Rhi- 
nolophus lepidus ( n = 36) and Myotis sp. (n = 1) were 
captured from nine locations throughout Singapore from 
April 2011 to March 2014. Collections were undertaken with 
approval from the National Parks Board (NP/RP11-011-3a) 
and the National University of Singapore Institutional Animal 
Care and Use Committee (IACUC Permit # B01/12). Bats 
were captured in-flight using a ground level or elevated mist 
net or a harp trap. Oral and rectal swabs were collected and 
treated in the same manner as the faecal and urine samples. 

Next-generation sequencing 

We performed next-generation sequencing on pooled faecal 
material and pooled urine collected from a colony of 
E. spelaea for virus discovery. Ten grams of faecal matter 
pooled from collections taken from plastic sheets laid under 
the colony on 14 March, 28 March and 11 April 2013 was 
vortexed in TBS buffer (25 him Tris, 150 him NaCl and 
Roche ULTRA protease inhibitor cocktail tablet) to resus¬ 
pend the material into solution. The samples were then 
homogenized using silica beads (MP FastPrep - 24 bead 
beater; MP Biomedicals, Santa Ana, CA, USA). After 
homogenization, the samples underwent a series of cen¬ 
trifugations where the supernatant was transferred to new 
tubes (full protocol available upon request). After the final 
centrifugation step, the supernatant was discarded and the 
pellet was resuspended in 500 /A of TBS buffer with pro¬ 
tease inhibitors and kept at — 80°C. A library was prepared 
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from the faecal processing and the serological capture 
extracted RNA/DNA. The SMARTer Universal Low Input 
RNA Kit (Clontech Laboratories, Mountain View, CA, 
USA) used with low concentration nucleic acids. The dou¬ 
ble-stranded cDNA library was purified with AMPure beads 
(Agencourt, Beckman Coulter, Brea, CA, USA), adaptors 
were removed with a digestion mix, and the library was 
amplified with an Advantage 2 PCR kit (Clontech Labora¬ 
tories). The libraries were analysed on a bioanalyser (Agi¬ 
lent Technologies, Santa Clara, CA, USA), and the reaction 
was run on both the Illumina HiSeq and MiSeq machines. 

Urine samples from an E. spelaea colony in Singapore on 
24 April, 8 May and 20 May 2014 were pooled and cen¬ 
trifuged at 10 000 g for 3 min. TRI-Reagent was added to 
urine viral supernatant and RNA extracted using Direct- 
zol™ RNA MiniPrep (#R2050; Zymo Research Corpora¬ 
tion, Irvine, CA, USA) as per manufacturer’s instruction. 
RNA was subjected to in-column DNase I (#M0303S; New 
England BioLabs Inc., Ipswich, MA, USA) digestion as 
instructed in Direct-zol™ RNA MiniPrep manual. 
Extracted and DNase I digested RNA was treated with 
Ribo-Zero™ Gold rRNA Removal Kit (Epidemiology) 
(#MRZE706; Epicentre, Madison, WI, USA) as described in 
the instruction manual. 

Two 500-bp next-generation sequencing cDNA libraries 
(urine-MiSeq-25 and urine-MiSeq-27) were constructed 
from the same sample using NEBNext® Ultra™ Directional 
RNA Library Prep Kit for Illumina® (#E7420S; New Eng¬ 
land BioLabs Inc.) as per manufacturer’s instruction. cDNA 
libraries were visualized on 1.5% agarose gel, excised and 
purified using Zymoclean Gel DNA Recovery Kit (#D4007; 
Zymo Research Corporation) as instructed. 

Fastq read files generated from the Illumina MiSeq and 
HiSeq runs were analysed using the DIAMOND alignment 
tool (Buchfink et al., 2015), against reference sequences in 
a custom protein database generated from RefSeq virus and 
neighbour nucleotide records (http://www.ncbi.nlm.nih.- 
gov/genome/viruses/) (Brister et al., 2015) using a custom 
Perl script. DIAMOND alignment results were imported 
into MEGAN 6 (Huson and Mitra, 2012), using a naive 
LCA algorithm to perform taxonomic assignment and visu¬ 
alization of aligned reads. All reads taxonomically assigned 
to family Coronaviridae were extracted and pooled, and 
SPAdes 3.6.2 (Bankevich et al., 2012) was used to perform 
de novo assembly of pooled reads. 

PCR screening 

Oral and rectal swabs or fresh faeces from each individual 
were pooled and screened for coronaviruses. A total of 974 
urine samples (947 samples from individuals and 27 sam¬ 
ples pooled by collection period) collected biweekly from 
the E. spelaea colony over 4 years from May 2011 to May 
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2015 in addition to 150 fresh faecal pellets from the above- 
mentioned bat colony collected during 3 time points over 
6 weeks (April-May 2014) were included. RNA was 
extracted from these samples using the QIAxtractor auto¬ 
mated purification system (QIAGEN, Hilden, Germany) 
following the manufacturer’s instructions. cDNA was 
synthesized using Superscript II Reverse Transcriptase 
(Invitrogen, Carlsbad, CA, USA) and screened using coron- 
avirus family-specific primer (PanCor IN-6 5'- GGTTGG 
GACTATCCTAAGTGTGA -3'; PanCor IN-7 5'- CCATCA 
TCAGATAGAATCATCATA -3') (Drosten et al., 2003) that 
targeted the RNA-dependent RNA polymerase gene. The 
PCR protocol used an initial denaturation step of 95°C for 
10 min, 40 cycles of denaturation at 95°C for 1 min, 
annealing at 56°C for 1 min, extension at 72°C for 1 min, 
a final extension at 72°C for 5 min. PCR products were 
screened on a 1.5% agarose gel and visualized on a gel doc 
machine. RNA from cultured Coronaviruses 229E, OC43 
and NL63 was extracted manually, and cDNA was synthe¬ 
sized to serve as positive controls. 

Phylogenetic analysis 

Coronavirus sequences representing diverse groups were 
downloaded from GenBank for three genes used in this 
study: RdRp (RNA-dependent RNA polymerase genes), Hel 
(helicase: nspl3), E (envelope) and N (nucleocapsid), 
including a-coronavirus, p-coronavirus (comprised of lin¬ 
eages A-D), and MERS-coronavirus, obtained from a wide 
range of hosts (bats, felines, human, avian and camels). Two 
Indonesian bat coronavirus that were recently detected from 
Moluccan naked-backed fruit bats (Dobsonia moluccensis) 
were also included (Stamatakis, 2014; Anindita et al., 2015). 
Partial sequences of the novel bat coronavirus generated 
were aligned with this data set and alignments manually 
edited in Geneious v9.0.3 (Biomatters, Auckland, New Zeal¬ 
and). Final aligned data sets comprise 82 RdRp sequences 
(2814 bp in length), 69 Hel sequences (1827 bp), 67 E 
sequences (240 bp) and 49 N sequences (578 bp) that were 
analysed separately (see accession numbers in Table SI). 
Phylogenetic trees were reconstructed using maximum likeli¬ 
hood implemented in RAxML v8.0.14 (Scientific Computing 
Group, Heidelberg Institute for Theoretical Studies, Heidel¬ 
berg, Germany) (Stamatakis, 2014), and branch support was 
assessed with 1000 nonparametric bootstrap replicates, with 
only values greater than 50% indicated at the major nodes. 
Trees were rooted with one representative of the gamma- 
coronaviruses (y-CoV, accession number: AY338732). 

Host distribution mapping 

Bat species distribution and country shapefiles were down¬ 
loaded from the International Union for Conservation of 
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Nature (IUCN, 2016) and imported into qGIS to create a 
map displaying the distribution of known Betacoronavirus 
lineage D reservoirs from South Asia and South-East Asia 
(QGIS Development Team, 2016). A map was also made of 
the overlapping distribution of the three cave roosting bat 
species. 

Results 

The Illumina NGS data sets from libraries generated from 
E. spelaea faeces and urine produced a total of 772 coron- 
avirus reads (Table 1). The libraries constructed from the 
urine samples produced 712 coronavirus reads, while the 
HiSeq and MiSeq runs from the faecal samples generated 60 
reads. The MiSeq run on the pooled faeces produced the 
fewest reads. The read lengths were from a minimum of 
76 bp to a maximum of 251 bp, with a mean of 240 bp 
(stdv 34.7 bp). A total of nine contigs were generated from 
the de novo assembly (min. 392 bp and max. 1092 bp). All 
431 pooled oral-rectal samples from the six bat species were 
PCR-negative. The urine and faecal samples (n = 1124) col¬ 
lected from the E. spelaea colony were also PCR-negative. 

We analysed four individual data sets ( RdRp , Hel , E and 
N) in the phylogenetic analysis, with each data set including 
alpha- and beta-coronaviruses. Maximum-likelihood phy- 
logenies show Betacoronavirus is composed of four distinct 
and well-supported lineages (A-D), although the intra¬ 
generic relationships of beta-coronaviruses vary between 
the four gene trees (Figs 1 and 2). The RdRp phylogeny 
places lineage C basal to lineages A, B and D, which form a 
monophyletic group with only moderate support (boot¬ 
strap, BS = 58%). This interlineage relationship is consis¬ 
tent with a recent bat CoV study from Thailand 
(Wacharapluesadee et al., 2015). In contrast, in the Hel 
phylogeny, lineage A is basal to the monophyletic lineages 
B, C and D (BS = 88%), with lineages B and D forming a 
strongly supported sister group (BS = 90%). This intra¬ 
generic relationship based on the Hel gene is also concor¬ 
dant with previous studies (Lau et al., 2012; Woo et al., 
2012b). The E gene phylogeny resembles that of the RdRp 
in grouping lineages A and D as sister taxa but with low 
support (BS = 51%). We also compared the genetic diver¬ 
gence of the RdRp gene within the different lineages: lineage 
D exhibits a higher level of sequence variation with only 
58.4% nucleotide identity, in contrast to lineages A 


Table 1 . Total number of coronavirus reads from four NGS data sets 
on pooled faeces and urine from a colony of Eonycteris spelaea 



Total reads 

Total corona reads 

Urine MiSeq 

18 547 895 

712 

Faeces HiSeq 

68 584 413 

26 

Faeces MiSeq 

4 952 973 

34 
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(83-100% nucleotide identity), B (92-100%) and C (82- 
100%) that show a higher percentage similarity. The N gene 
phylogeny indicated lineage D is basal to other lineages A, 
B and C; however, this interlineage relationship was statisti¬ 
cally insignificant. 

The RdRp , Hel , E and N gene phylogenies all place the 
novel CoV identified from E. spelaea unequivocally in the 
well-supported monophyletic Betacoronavirus lineage D 
(RdRp and Hel : 100% BS; E: 91% BS; N: 97%). The HKU9 
prototype virus for lineage D Betacoronavirus was reported 
from Rousettus leschenaultii in Hong Kong in 2005 (Woo 
et al., 2007). Within lineage D, partial RdRp sequence of 
the novel CoV possesses 77.2-97.4% nucleotide sequence 
identity with the remaining lineage D bat-CoVs. In con¬ 
trast, the partial Hel and E genes of the E. spelaea CoV 
sequences exhibit greater similarity (Hel: 80-82.5%; E: 
70.6-97.9%) with other lineage D viruses, although this 
may due to the limited number of lineage D Hel and E gene 
sequences that are available. The N gene indicated the CoV 
of E. spelaea showed more distant relationship with the 
remaining lineage D bat-CoVs (52-60% similarity). 

The RdRp phylogeny indicates that the E. spelaea CoV is 
most closely related to a virus from a large Asian roundleaf 
bat (Hipposideros lekaguli ), strongly supported with a boot¬ 
strap value of 100%. These two sequences form a sister 
group (bootstrap 93%) with coronaviruses from Cynop- 
terus sphinx , D. moluccensis , Hipposideros commersoni , Pte- 
nochirus jagori , R. leschenaultii and Rousettus aegypticus 
that comprise the majority of lineage D betacoronavirus 
sequences available. 

Discussion 

This is the first coronavirus detected from a bat in Singa¬ 
pore, and our study provides the first evidence of cave nec¬ 
tar bat E. spelaea harbouring a lineage D betacoronavirus. 
Coronaviruses were previously detected in E. spelaea in the 
Philippines, but no sequence data were generated (Watan- 
abe et al., 2010). Eonycteris spelaea has also previously been 
reported to harbour a paramyxovirus in China (Yuan et al., 
2014), Phnom-Penh bat virus (Queen et al., 2015), Issyk- 
kul virus (Calisher et al., 2006) and were also seropositive 
for Nipah virus in Malaysia (Yob et al., 2001). 

In South-East Asia, lineage D betacoronaviruses (also 
referred to as Ro-BatCoV HKU9: Woo et al., 2012a,b) have 
been detected in eight different bat species from three bat 
families belonging to two suborders, Yinpterochiroptera 
and Yangochiroptera: C. sphinx , D. moluccensis , E. spelaea , 
H. lekaguli , P. jagori , R. leschenaultii, Scotophilus heathi 
and Scotophilus kuhli (Woo et al., 2007; Tsuda et al., 2012; 
Anindita et al., 2015). A next-generation sequencing 
approach in Yunnan province China detected HKU9 coro¬ 
navirus sequences in a community dominated by 
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NC009657|Scofop/i//us /tu/]//7|Hainan|BtCoV/512|2005 
■ NC003436|Porcine|PEDV CV777 

EF203066|R/i/no/op/ius s/n/cus|HongKong|Bat-CoV-HKU2/HK/298|2004 
EF20m7\Rhinolophus s/n/cus|HongKong|Bat-CoV-HKU2/HK/33|2004 
NC002645|Homo sap/ens|HCoV-229E inf-1 
DQ4459111 Homo sap/'ens|Netherlands|HCoV-NL63|Amsterdam057|2004 
DQ445912|Homo sapiens\ Netherlands!HCoV-NL63|Amsterdam496|2003 
AY518894|Homo sap/'ens|Netherlands|HCoV-NL63 NL 
DQ848678|Feline|United Kingdom|FIPV ClJe 
100| DQ286389|Feline|FIPV DF-2 

NC007025|Feline|USA|FIPV WSU-79/1146 
AY994055|Feline|USA|FIPV 79-1146 
DQ811787|Porcine|USA|PRCV ISU-1 
DQ201447|Porcine|TGEV TS 
DQ811786|Porcine|USA|TGEV MillerM60| 1987 
DQ443743|Porcine|TGEV SC-Y 
DQ811789|Porcine|USA|TGEV virulentPurdue 
AJ271965|Porcine|USA|TGEV PUR46-MAD 
DQ811788|Porcine|TGEV Purdue-P115 
NC008315|7y/onycferis pachypus[Guangdong|Bat-CoV-HKU4/133|2005 
EF065508 |Tylonycteris pac/jypus|Guangdong|Bat-CoV-HKU4-4/LMH1f|2005 
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100P EF065509|Ptp/sfre//us afcramas|Guangdong|Bat-CoV-HKU5-1/LMH3f|2004 
EF06551 HPipistrellus af>ramus|Guangdong|Bat-CoV-HKU5/TT06f|2004 
EF065510|P/p/sfre//as a6ramus|Guangdong|Bat-CoV-HKU5-2/TT03f|2004 
EF065512|P/p/sfre//us abramus|Guangdong|Bat-CoV-HKU5/TT07f|2004 
KC776174[Homo sap/ens|Jordan|betacoronavirus 2c Jordan-N3|2012 
KM027258 |Homo sap/ens|Saudi Arabia|MERS-Jeddah-C8826/KSA|2014 > 

KR011263|Homo sap/ens|Saudi Arabia|MERS-Hu/Riyadh-KSA-2345|2015 o 

KR011265|77omo sap/'ens|Saudi Arabia|MERS-Hu/Riyadh-KSA-2466|2015 O 

KJ156952|Homo sap/ens|Saudi Arabia|MERS-Riyadh-4|2013 ^ 

KF186564|Homo sap/ens|Saudi Arabia|MERS-AI-Hasa-4|2013 
KF958702 |Homo sap/ensjsaudi Arabia|MERS-Jeddah-human-1|2013 LU 

KC667074 |Homo sap/ens|United Kingdom|betacoronavirus 2c England/Qatar|2012 
KC164505|Homo sap/'ensjUnited Kingdom|betacoronavirus England 112012 
NC019843|Homo sap/'ens|Netherlands|MERS-HCoV-EMC|2012 
KP719932|Came/us dromecfar/uslUnited Arab Emirates|Camel/UAE/D1243|2014 
DQ648856 |Rhinolophus fe/rameg(j/™m|Hubei|SARS-bat-BtCoV/273|2004 
NC009695 |Rhinolophus fe™meg™i/m|Hubei|SARS-bat-Rf112004 
HC00%98\Rhinolophus pearson/|Guangxi|SARS-bat-Rp3|2004 
NC009696|Rt)/'no/opt)us macrof/s|Hubei|SARS-bat-Rm1|2004 
DQ648857|Rh/>70/op/ius macro//s|Hubei|SARS-bat-BtCoV/279|2004 
DQ084199|Rb/noiop/)t/s s/n/cus|Hong Kong|SARS-bat-HKU3-2|2005 
DQ084200|R/r/r)o/op/]us s/m'cus|Hong Kong|SARS-bat-HKU3-3|2005 
NC00A718\Rhinolophus s/n/cus|Canada|SARS-bat-HKU3-1/Tor2|2003 
100 I KJ020604|Cynopferus sp/)/77x|Thailand|BtCoV/B55700-3/Cyn sph/CB/Tha/5|2012 
KJ020607 |Scotophilus f)eaff)//|Thailand|BtCoV/B55762/S hea/CB/Tha/6|2012 
1001 Kmom\Scotophilus /a/W//|Thailand|BtCoV/B55740/S kuh/CB/Tha/6[2012 
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EF065514|Rouseffi/s /esct)enaa/f/7|Guangdong|Bat-CoV-HKU9-2 BF 017I|2005 
100r AB918718|Dobsom'a mo/uccens/s|Gorontalo|Bat-CoV-IFB2012-17F|2012 
AB918719|Dobson/a mo/uccens/s|Gorontalo|Bat-CoV-IFB2012-8F|2012 
EF065515|Rousef/us /esc/?enau///7|Guangdong|Bat-CoV-HKU9-3 BF 493112005 
81 KJ020605|Cynopferus spt)/ 77 x|Thailand|BtCoV/B 55735 /C sph/CB/Tha/5|2012 


100 



100 


99 


91 


0.3 


100 


100 



Coronavirinae 

■ a-CoV 

(3-CoV lineage A 
p-CoV lineage B 
(3-CoV lineage C 

■ p-CoV lineage D 


AB683970|Pfenoct)/7as y'agror/'| Philippines|Bat-CoV-2265|2010 
GQ920804|Rouseffus aegypt/acus|Kenya|Bat-CoV-BtKY05|2006 
■ GQ920806 |Hipposideros commer'sof7/|Kenya|Bat-CoV-BktY07|2006 
I GQ920825|Rousettus aegypf/acus|Kenya|Bat-CoV-BtKY28|2006 
8 iW EF065516|Rouseftus/esc/)enau/t/7|Guangdong|Bat-CoV-HKU9-4 BF 141112005 
54 L EF065513|Rouseffus/esct)enaa/f//|Guangdong|Bat-CoV-HKU9-1 BF 005112005 
DQ415906 |Homo sap/ens|Hong Kong|CoV-HKU1A-N9[2003 
DQ415901 |Homo sap/ens|Hong Kong|CoV-HKU1A-N24|2003 
DQ415905|/7omo sap/ens|Hong Kong|CoV-HKU1A-N7|2003 
DQ415914|Homo sap/'ens|Hong Kong|CoV-HKU1A-N18|2003 
DQ415908|Homo sap/ens|Hong Kong|CoV-HKU1 A-N1112003 
NC006577[Homo sap/ens|Hong Kong|CoV-HKU1A-N1|2003 
DQ415903 |Homo sap/ens|Hong Kong|CoV-HKU1A-N3|2003 
AY903460|Homo sap/'ens|Belgium|HCoV-OC43|2004 
NC007732|Porcine|Belgium|PHEV VW572 
U00735|Bovine|BCoV Mebus 
AF220295|Bovine|BCoV Quebec 
DQ811784|Bovine|BCoV DB2|1983 
EF424620|Bovine|USA|BCoV R-AH187 
DQ915164|Alpaca|USA|BCoV NA 
NC003045|Bovine|BCoV-ENT 
EF424615]Bovine|USA|BCoV E-AH65 
AF391542|Bovine|BCoV-LUN 
EF424618|Bovine|USA|BCoV R-AH65-TC 
EF424617|Bovine|USA|BCoV R-AH65 



100 


100 


EF203067 Bat-CoV-HKU2 HKU2/HK/33/2004 
EF203066 Bat-CoV-HKU2 HKU2/HK/298/2004 
NC003436|Porcine|PEDV CV777 
NC009657|Scofopt)/7us /tuW/7|Hainan|BtCoV/512|2005 
NC002645|Homo sap/ens|HCoV-229E inf-1 
AY518894|Homo sap/ens|Netherlands|HCoV-NL63 NL 
DQ445912|Homo sap/ens|Netherlands|HCoV-NL63|Amsterdam496|2003 
DQ445911 \Homo sap/ens|Netherlands|HCoV-NL63|Amsterdam057|2004 
DQ848678|Feline|United Kingdom|FiPV ClJe 
AY994055|Feline|USA|FIPV 79-1146 
DQ286389|Feline|FIPV DF-2 
NC007025|Feline|USA|FIPV WSU-79/1146 
DQ811787|Porcine|USA|PRCV ISU-1 
DQ811786|Porcine|USA|TGEV MillerM60| 1987 
DQ201447|Porcine|TGEV TS 
DQ443743|Porcine|TGEV SC-Y 
DQ811788|Porcine|TGEV Purdue-P115 
AJ271965|Porcine|USA|TGEV PUR46-MAD 
DQ811789|Porcine|USA|TGEV virulentPurdue 
■ AY884001|Homo sap/ens|Hong Kong|CoV-HKU1B N2|2003 
NC006577|Homo sap/ens|Hong Kong|CoV-HKU1A-N1|2003 
DQ415904|Homo sap/ens|Hong Kong|CoV-HKU1A-N6|2003 
DQ4159011 Homo sap/er?s|Hong Kong|CoV-HKU1A-N24|2003 
DQ415896|Homo sap/ens|Hong Kong|CoV-HKU1A-N19|2003 
DQ415900|Homo sap/ens|Hong Kong|CoV-HKU1A-N23|2003 
DQ415914|Homo sap/ens|Hong Kong|CoV-HKU1A-N18|2003 
NC007732|Porcine|Belgium|PHEV VW572 
U00735|Bovine|BCoV Mebus 
AF220295|Bovine|BCoV Quebec 
DQ811784|Bovine|BCoV DB2|1983 
AF391542|Bovine|BCoV-LUN 
DQ915164|Alpaca|USA|BCoV NA 
NC003045|Bovine|BCoV-ENT 
EF424620|Bovine|USA|BCoV R-AH187 
KX452688|£onycferis spe/aea|Singapore|Bat-CoV|2014 ★ 

AB918719 |Dobsonia mo/uccef 7 s/s|Gorontalo|Bat-CoV-IFB 201 2-8F|2012 
AB918718|Dobson/'a mo/tvcce/7s/'s|Gorontalo|Bat-CoV-IFB2012-17F|2012 

EF065515|Rouseffus /esc/)enau/t/7|Guangdong|Bat-CoV-HKU9-3 BF 493112005 
EF065514jRouseffus /esct)enau/t//|Guangdong|Bat-CoV-HKU9-2 BF 017I|2005 
EF065516|Rouseftus /esct)er 7 au/t/ 7 |Guangdong|Bat-CoV-HKU 9-4 BF 141112005 
69*- EF065513|Rousef/us /esc/ienau///7|Guangdong|Bat-CoV-HKU9-1 BF 005112005 
UC00%$3\Rhinolophus pearson/|Guangxi|SARS-bat-Rp3|2004 
100 » DQ084199|R/r/no/op/)us s/f 7 /cus|Hong Kong|SARS-bat-HKU3-2|2005 
DQ084200\Rhinolophus s/n/cus|Hong Kong|SARS-bat-HKU3-3|2005 
NC004718|RWno/op/)us s/mcus|Canada|SARS-bat-HKU3-1/Tor2|2003 
NC009695|R/i/no/op/)us fe/rumegu/™m|Hubei|SARS-bat-Rf1|2004 
•gDQ648856|Rb(7]o/op/ius fe™megu/7ium|Hubei|SARS-bat-BtCoV/273|2004 
85|j DQM8857\Rhinolophus macrof/'s|Hubei|SARS-bat-BtCoV/279|2004 
100* UC0096%\Rhinolophus macrof/s|Hubei|SARS-bat-Rm1|2004 
EF065509|P/p/sfre//us abramus|Guangdong|Bat-CoV-HKU5-1/LMH03f|2004 
EF065510|P/p/'sfre//L/s abramus|Guangdong|Bat-CoV-HKU5-2//TT03f|2004 
EF065511 |P/p/sfre//us at>ramus|Guangdong|Bat-CoV-HKU5/TT06f|2004 
EF065512|P/p/sfre//us abramus|Guangdong|Bat-CoV-HKU5/TT07f|2004 
NC008315| Tylonycteris pac/jypus|Guangdong|Bat-CoV-HKU4/133|2005 
NC009019|Ty/onycteris pac/)ypus|Guangdong|Bat-CoV-HKU4-1/B04f|2005 
EF065507| Tylonycteris pac/jypus|Guangdong|Bat-CoV-HKU4-3/B07f|2005 
EF065508| Tylonycteris pac/jypus|Guangdong|Bat-CoV-HKU4-4/LMH1f|2005 
EF065506| Tylonycteri spac/rypus|Guangdong|Bat-CoV-HKU4-2/B05f|2005 
KM027258 |Homo sapiens |Saudi Arabia|MERS-Jeddah-C8826/KSA|2014 
NC019843|Homo sap/ens|Netherlands|MERS-HCoV-EMC|2012 
KC776174|Homo sap/'ens|Jordan|betacoronavirus 2c Jordan-N3|2012 
KR011265 |Homo sap/ens|Saudi Arabia|MERS-Hu/Riyadh-KSA-2466|2015 
KR011263|Homo sapiens |Saudi Arabia|MERS-Hu/Riyadh-KSA-2345|2015 
KP719932|Came/us dromeda™s|United Arab Emirates|Camel/UAE/D1243|2014 
KF186564|Homo sap/ens|Saudi Arabia|MERS-AI-Hasa-4|2013 
KF958702 |Homo sap/ens|Saudi Arabia|MERS-Jeddah-human-1|2013 
KC164505|Homo sap/ens|United Kingdom|betacoronavirus England 112012 
KJ156952|Homo sap/'ens|Saudi Arabia|MERS-Riyadh-4|2013 
KC667074 |Homo sap/'ens|United Kingdom|betacoronavirus 2c England/Qatar|2012 
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Fig. 1. Phylogenetic trees of the RdRp and Hel nucleotide sequences of alpha- and beta-coronaviruses reconstructed using maximum-likelihood 
method in RAxML, with 1000 bootstrap replicates. Bootstrap support values greater than 50% are indicated at major nodes. The trees were rooted 
with one representative of the gammacoronavirus (y-CoV, accession number: AY338732). Coloured branches denote different betacoronavirus lin¬ 
eages (i.e. A-D). New CoV sequence collected from a bat species (Eonycteris spelaea) in this study is marked by an asterisk. The scale bar indicates the 
number of nucleotide substitutions per site. 


Hipposideros armiger , but the host identity of the faeces 
tested was not confirmed (Ge et al., 2012). The bat species 
known to harbour lineage D Betacoronavirus are widely 
distributed across South Asia and South-East Asia, and 
there are large areas where these distributions overlap 
(Fig 3, Table S2). Lineage D betacoronaviruses were also 
detected in two bat species in Kenya, Rousettus aegyptiacus 
and H. commersoni , and one in Madagascar, Pteropus rufus 
(Tong et al., 2009; Razanajatovo et al., 2015). 

Interestingly, the pooled urine samples where coron- 
avirus was detected in the Illumina NGS were negative 
with traditional RT-PCR, indicating a low amount of 
virus. These screening results possibly indicate that there 
is a low prevalence in the bats in Singapore, low virus 


titres in the samples and/or we missed shedding foci dur¬ 
ing our sampling. The only bat species we detected coro- 
naviruses in were the cave nectar bats (E. spelaea ), and 
this was using deep sequencing where the majority of the 
reads were detected in the urine sample. We did not 
detect this in the other insectivorous or frugivorous bat 
species with this well-validated primer set; however, we 
did not perform NGS on samples from other species in 
Singapore. The majority of coronavirus detections have 
been from insectivorous bats, but this may have been 
spurred by the discovery of SARS-like CoV in Rhinolo¬ 
phus bat species in China (Drexler et al., 2014). Different 
coronaviruses can infect the same species of bats, even 
very divergent viruses (Wacharapluesadee et al., 2015). 
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NC009657|Scofop/i//us toW//|Hainan|BtCoV/512|2005 
NC003436|Porcine|PEDV CV777 
NC002645|Homo sap/ens|HCoV-229E inf-1 

100 | EF203067|R/)/no/ophus sm/cus|HongKong|Bat-CoV-HKU2/HK/33|2004 
EF203066|R/)/no/opfws sm/cus|HongKong|Bat-CoV-HKU2/HK/298|2004 
10 o| AY518894|Homo sap/ens|Netherlands|HCoV-NL63 NL 

DQ445911 |Homo sap/ens|Netherlands|HCoV-NL63|Amsterdam057|2004 
DQ445912|Homo sap/e/7s|Netherlands|HCoV-NL63|Amsterdam496|2003 
DQ286389|Feline|FIPV DF-2 
NC007025|Feline|USA|FIPV WSU-79/1146 
AY994055| Feline|USA| FIPV 79-1146 
DQ848678|Feline|United Kingdom|FIPV ClJe 
DQ811788|Porcine|TGEV Purdue-P115 
DQ443743|Porcine|TGEV SC-Y 
AJ271965|Porcine|USA|TGEV PUR46-MAD 
DQ811789|Porcine|USA|TGEV virulentPurdue 
DQ811787|Porcine|USA|PRCV ISU-1 
DQ811786|Porcine|USA|TGEV MillerM60|1987 
DQ201447 |Porcine|TGEV TS 

EF065514|Rouseffus /eschenau/(/7|Guangdong|Bat-CoV-HKU9-2 BF 017I|2005 
KX452686|Eonycfer/s spe/aea|Singapore|Bat-CoV|2014* 

EF065513|Rouseffus /eschenau/ft7|Guangdong|Bat-CoV-HKU9-1 BF 005112005 | 
p EF065515|Rousefft/s /esc/ienau/f//|Guangdong|Bat-CoV-HKU9-3 BF 493l|2005 
L EF065516|Rouseffus /eschenau/f/7|Guangdong|Bat-CoV-HKU9-4 BF 1411|2005 
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Coronavirinae 

■ a-CoV 

(3-CoV lineage A 
p-CoV lineage B 
|3-CoV lineage C 
P-CoV lineage D 


NC006577|Homo sap/ens|Hong Kong|CoV-HKU1 A-N112003 
DQ415903|Homo sap/e/?s|Hong Kong|CoV-HKU1 A-N3|2003 
DQ415905|Homo sap/ens|Hong Kong|CoV-HKU1 A-N7|2003 
DQ415908|Homo sap/ens|Hong Kong|CoV-HKU1 A-N1112003 
DQ415901|Homo sap/e/is|Hong Kong|CoV-HKU1A-N24|2003 
DQ415906|Homo sap/e/?s|Hong Kong|CoV-HKU1 A-N9|2003 
DQ415914|Homo sap/e/?s|Hong Kong|CoV-HKU1 A-N18|2003 
NC007732|Porcine|Belgium|PHEV VW572 
AF391542|Bovine|BCoV-LUN 
EF424620|Bovine|USA|BCoV R-AH187 
NC003045|Bovine|BCoV-ENT 
DQ915164|Alpaca|USA|BCoV NA 
U00735|Bovine|BCoV Mebus 
AF220295|Bovine|BCoV Quebec 
DQ811784|Bovine|BCoV DB2|1983 
NC009695|R/j/7jo/opPus ferrumegui7]t/m|Hubei|SARS-bat-Rf1|2004 
DQ648856| Rhinolophus fe™megu/nt/m|Hubei|SARS-bat-BtCoV/273|2004 
NC009693|Rhi7io/opPi/s peareoni|Guangxi|SARS-bat-Rp3|2004 
NC00%%\Rhinolophus macrof/s|Hubei|SARS-bat-Rm1|2004 
DQ648857|R/j/no/op/jus macrof/s|Hubei|SARS-bat-BtCoV/279|2004 
DQ084199|Rfimo/op/jus s/mcus|Hong Kong|SARS-bat-HKU3-2|2005 
DQ084200|Rfiino/op/jt/s s/n/eus|Hong Kong|SARS-bat-HKU3-3|2005 
NC004718|Rfi/7jo/opfius s/m'cus|Canada|SARS-bat-HKU3-1/Tor2|2003 
EF0655091P/p/sfre//[/s abra/nus|Guangdong|Bat-CoV-HKU5-1/LMH03f|2004 
EF065511 \Pipistrellus aPramus|Guangdong|Bat-CoV-HKU5/TT06f|2004 
EF0655101 Pipistrellus abramus|Guangdong|Bat-CoV-HKU5-2/TT03f|2004 
EF065512|P/p/sfre//us abra/nus|Guangdong|Bat-CoV-HKU5/TT07f|2004 

EF065508| Tylonycteris pacbypus|Guangdong|Bat-CoV-HKU4-4/LMH1f|2005 
NC008315|7y/onycferis pacbypus|Guangdong|Bat-CoV-HKU4/133|2005 
EF065506|7y/onycferis pacbypus|Guangdong|Bat-CoV-HKU4-2/B05f|2005 
NC009019| Tylonycteris pacftypus|Guangdong|Bat-CoV-HKU4-1/B04f|2005 
EF065507|7y/onycferis pacftypus|Guangdong|Bat-CoV-HKU4-3/B07f|2005 
KR011265|Homo sap/ens|Saudi Arabia|MERS-Hu/Riyadh-KSA-2466|2015 
KC164505|Homo sap/ens|United Kingdomjbetacoronavirus England 1|2012 
KC776174|Homo sap/ens|Jordan|betaeoronavirus 2c Jordan-N3|2012 
KJ156952|Homo sap/ens|Saudi Arabia|MERS-Riyadh-4|2013 
KP719932|Came/us tfromec(arius|United Arab EmiratesjCamel/UAE/D1243|2014 
KR011263|Homo sap/ens|Saudi Arabia|MERS-Hu/Riyadh-KSA-2345|2015 
KC667074| Homo sap/ens|United Kingdom|betacoronavirus 2c England/Qatar|2012 
KF958702|Homo sap/e/7s|Saudi Arabia|MERS-Jeddah-human-1|2013 
KM027258I Homo sapiens |Saudi Arabia|MERS-Jeddah-C8826/KSA|2014 
KF186564jHomo sap/ens|Saudi Arabia|MERS-AI-Hasa-4|2013 
NC019843|Homo sap/ens|Netherlands|MERS-HCoV-EMCi2012 


100 


JN825712|Porcine| PEDV-B J-2011-112011 
FJ710055| Hipposideros sp. |Ghana|Bat-CoV/344|2008 
DQ243942|Homo sap/ens|Australia|HCoV-229E-16/6/82 
KT368911 |Camel|Saudi Arabia|Camel/Jd8|2015 
JQ410000|Alpaca|USA|Alpaca/CA-08-0112008 
FJ710056|H/ppos/deros sp. |Ghana|Bat-CoV/19|2008 
100 I DQ445911|Homo sapiens|Netherlands|HCoV-NL63/Amsterdam057|2004 
l DQ462768|Homo sapiens|Netherlands|HCoV-NL63/246|2003 
— DQ666344|IW/n/opferus sp.|Hong Kong|Bat-CoV-HKU8/AFCD77|2006 
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' EF203064I Rhinolophus sp.|Hong Kong|Bat-CoV-HKU2/GD/430|2006 
■ HQ728486|C/iaerephon sp.|Kenya|Bat-CoV/KY22|2006 
■ JQ989271 IRouseffus sp. |China| Bat-CoV-H KU10/175A|2005 
99r FJ917524|Feline|USA| FI PV UCD14|2007 
100 I — GU017120|Feline|ltaly|FIPV Ml/cat8/04|2004 
■ KP202848|Porcine|China|TGEV SHXB|2013 
• AF104420|Porcine|United Kingdom|TGEV 96-1933 

KX452687|Eonycteris spe/aea|Singapore|Bat-CoV|2014* 

HQ728482|E/d0/0f? be/i/am|Kenya|Bat-CoV-BtKY24|2006 
1 AB543561|Cynopferus brac/iyof/s|Philippines|Bat-CoV-Diliman1525G2|2008 
HM211101|Rousetfus /escdenau/t/7|Guangdong|Bat-CoV-HKU9-10-2 BF 506l|2006 
EF065516|Rouseffus /esc/jenau/f//|Guangdong|Bat-CoV-HKU9-4 BF 1411|2005 
EF065515|Rouseffus /esc/]enau/f//|Guangdong|Bat-CoV-HKU9-3 BF 493l|2005 
10o| HM211098|Rouse(fus/escdenau/f/7|Guangdong|Bat-CoV-HKU9-5-1 BF 258l|2005 
HM211100|Rousef(us /escPenau/f/7|Guangdong|Bat-CoV-HKU9-10-1|2006 
1001 HQ728483|Rouseftus aegypf/acus|Kenya|Bat-CoV-BtKY06|2006 

EF065513|Rouseftus /esc/ienau/f;/|Guangdong|Bat-CoV-HKU9-1 BF 005l|2005 
1001 KF686341|Homo spaiens|USA|HCoV-HKU1-10|2010 
” HM034837|Homo spa/ens|France|HCoV-HKU1/Caen1|2005 
AF058943|Bovine|USA|BCoV-LSU-94LSS-051-2 
KF572710|Homo spa/ens|China|FICoV-OC43/10241-10|2010 
KF923898|Homo spa/ens|China|HCoV-OC43/3184A|2012 
JN874562|Rabbit|China|HKU14-10|2007 
AF061835|Murine|MHV-2 

KF850449|Raffus sp. |USA|MHV/MG6644-Holmes10| 1996 

JQ989269|Hipposideros sp.|Hong Kong|SARS-bat-CoV-HKU10|2005 
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GQ153545\Rhinolophus sp. |China|SARS-bat-CoV-HKU3-10|2006 
DQ071614|Rto'no/opf)us sp.|China|SARS-bat-Rp2|2004 
_ ■ FJ588688|Rbmo/op/]us s/nicus|China|SARS-Bat-CoV-Rs806|2006 
EF199653|China|SARS-TJ01 
AY283795|Homo sapiens|Singapore|SARS-Sin2677 
AY283798| Homo sapiens|Singapore|SARS-Sin2774 
■ NC014470|Rfnno/ophus Was/7|Bulgaria|SARS-bat-CoV-BM48-31/BGR|2008 
EF065509|P/pisfre//us abramus|GuangdongjBat-CoV-HKU5-1/LMH03f|2004 
84 I » EF065510|P/p/sfre//us abramus|Guangdong|Bat-CoV-HKU5-2/TT03f|2004 
| 1001 EF065508| Tylonycteris pachypus|Guangdong|Bat-CoV-HKU4-4/LMH1f|2005 
' EF065506|7y/onycferis pacbypi/s|Guangdong|Bat-CoV-HKU4-2/B05f|2005 
KJ741003|Environment|Saudi Arabia|MERS-Jeddah-air-1|2014 
KT368871 |Camel|Saudi Arabia|MERS-Ridyadh/Ryl 62|2015 
KM609328|Camel|United Arab Emerites|MERS-AD-UAE-11146(20141 
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Fig. 2. Phylogenetic tree of the E and N nucleotide sequences of alpha- and beta-coronaviruses reconstructed using maximum-likelihood method in 
RAxML, with 1000 bootstrap replicates. Bootstrap support values greater than 50% are indicated at major nodes. The tree was rooted with one repre¬ 
sentative of the gammacoronavirus (y-CoV, accession number: AY338732). Coloured branches denote different betacoronavirus lineages (i.e. A-D). 
New CoV sequence collected from a bat species (Eonycteris spelaea) in this study is marked by an asterisk. The scale bar indicates the number of 
nucleotide substitutions per site. 


The close evolutionary relationship of lineage D betacoro¬ 
naviruses detected from geographically separated bat spe¬ 
cies indicates that the spread of lineage D infections in 
bats is not host-restricted. 

In Germany, bat parturition peaks appear correlated with 
bat coronavirus shedding (Drexler et al., 2011); it would be 
beneficial to determine the pregnancy cycles of bats in the 
tropics, which is largely unknown or based on studies with 
limited geographical scope. In the tropics, the constant 
availability of resources means there is less dependence on 
temporal windows of food. In Malaysia, E. spelaea has been 
found with dependent young in 8 of 12 months, indicating 
that individuals are reproducing throughout the year (King¬ 
ston et al., 2006), with two peak periods of pregnancy in 
June and September. In the Philippines, it was also found 
that E. spelaea has two annual peaks in reproduction 
although with some variation (Krutzsch, 2005). The first 


peak centred on March or April while the second centred 
on or near August (Heideman and Utzurrum, 2003). These 
bats typically roost in large numbers in caves, although the 
studied colony in Singapore roosts under an overpass and 
numbers over 3000 individuals (Lee, unpublished data). 
This opens up the possibility that there is a persistent source 
of virus in the colonies and immunologically naive individ¬ 
uals have the opportunity for exposure from infected bats 
as they are tightly packed and may fight for specific sites in 
a roosting area. This is especially true for larger and older 
males because it is hypothesized that E. spelaea exhibits a 
resource defence polygynous mating system, where males 
spend significantly more time and energy for roost defence 
and surveillance (Bumungsri et al., 2013). 

In addition, E. spelaea , R. leschenaultii and H. armiger 
co-roost in caves in India and in SE Asia and multispecies 
roosts may facilitate the transmission and recombination of 
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Diversity of Bat Betacoronaviruses in South-East Asia 



Fig. 3. Bat distribution map for each species in South-East Asia from which lineage D betacoronaviruses have been detected. These maps were gen¬ 
erated from IUCN maps in QGIS. 
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Table 2. Roosting and foraging behaviour of lineage D Betacoronavirus-positive bat species from South Asia and South-East Asia 


Bat species 

Roost site (cave/tree) 

Food (no. of 
plant species) 

Numbers in single roost 

References 

Eonycteris spelaea 

Caves in forested areas 

and man-made 

structure in urban areas 

Flowers/nectar 

Several thousands. 

Co-habits with 
other bat species 

(Francis et al., 2008) 

Cynopterus sphinx 

Foliage (tree and palms) only 

Fruits (wild and 
cultivated) 

3-7 bats per roost 

(Bates et al., 2008b) 

Dobsonia 

moluccensis 

Caves, sinkholes, boulder piles, 
old mines, disused buildings, 
dense vegetation 

Fruits 

Several thousands 

(Hutson et al., 2008) 

Hipposideros 

lekaguli 

Mainly limestone caves 

Insects 

300 individuals (per colony) 

(Bumrungsri per comm.) 

(Csorba et al., 2008) 

Rousettus 

leschenaultii 

Caves, old buildings, forts 
and disused tunnels 

Fruits and 

flowers/nectar 

A few to several thousand individuals 

(Bates and Helgen, 2008) 

Hipposideros 

armiger 

Caves and man-made structures 

Insects 

Hundreds of individuals. Co-habits 
with Rhinolophus and other bat species 

(Bates et al., 2008a) 

Ptenochirus 

jagori 

Mainly tree cavities, also caves and 
sheltered rock crevices 

Fruits (145), 
flowers/nectar 
(2) and leaves (7) 

Singly or in small groups in caves 

(Mickleburgh, 1992; 
Reiter and Curio, 2001; 
Ong et al., 2008) 

Scotophilus 

kuhli 

Tree cavities, foliage and man-made 
structures (in forest and 
human-modified environments) 

Insects 

(hymenopterans 
and dipterans) 

Several hundred individuals 

(Bates et al., 2008d) 

Scotophilus 

heathii 

Tree cavities, foliage and man-made 
structures (in forest and 
human-modified environments) 

Insects 

Singly or in colonies of up to 50 bats 

(Bates et al., 2008c) 


coronaviruses (Mendenhall, unpublished data; Struebig 
et al., 2005; Furey et al., 2011). Hipposideros lekaguli is 
thought to be associated with Rousettus species and even 
E. spelaea in SE Asia (Alviola et al., 2015). These species also 
have a large, overlapping distribution across SE Asia, which 
may facilitate the emergence of novel coronaviruses (Fig. 2). 
The foraging ecology of these bats E. spelaea is known to 
feed on a wide range of plant species and their high 
visitation frequency to the inflorescence of a few keystone 
plant species such as Durio zibethinus and Parkia speciosa 
(Heideman and Utzurrum, 2003; Krutzsch, 2005). Both 
C. sphinx and R. leschnaulti have also been found to feed on 
the same tree (Singaravelan and Marimuthu, 2004, 2008). 

The discovery of this novel coronavirus demonstrates the 
unsampled diversity of this virus family in bats. Betacoron¬ 
aviruses are responsible for two of the major coronavirus 
cross-species spillover events (SARS and MERS virus) from 
bats to incidental hosts, ultimately resulting in sustained 
transmission chains. The bat species where lineage D beta¬ 
coronaviruses have been detected share roosting sites, do so 
in large numbers and high density and also forage on the 
same plants. These may be factors that ultimately lead to the 
emergence of novel virus strains (Table 2). Even though the 
lineage D betacoronaviruses are not implicated as viruses 
capable of infecting humans, understanding their evolution 
can help us understand the evolution and ecological aspects 
of this medically important virus family. Furthermore, as 
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little is known about the evolutionary factors affecting the 
transmission of CoVs between hosts, additional studies, 
such as host cell receptor usage, are required to assess the 
potential risks of zoonotic transmission to humans. 
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